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The anisotropic superconductivity and the density wave have been investigated by apply- 
ing the Kadanoff- Wilson renormalization group technique to the quasi-one-dimensional system 
with finite-range interactions. It is found that a temperature (T) dependence of response func- 
tions is proportional to exp(l/r) in a wide region of temperature even within the one-loop 
approximation. Transition temperatures are calculated to obtain the phase diagram of the 
quasi-one-dimensional system, which is compared with that of the pure-one-dimensional sys- 
tem. Next-nearest neighbor interactions {V2) induce large charge fluctuations, which suppress 
the d^2_j^2-wave singlet superconducting (dSS) state and enhance the /-wave triplet supercon- 
ducting (/TS) state. From this effect, the transition temperature of /TS becomes comparable 
to that of dSS for large V2, so that fleld-induced /-wave triplet pairing could be possible. These 
features are discussed to comprehend the experiments on the (TMTSF)2PF6 salt. 
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1. Introduction 

The mechanism of superconductivity in quasi-one- 
dimensional (qld) systems has been attracting renewed 
interest since the appearance of the triplet supercon- 
ductivity suggested in the experiments on the qld or- 
ganic conductor (TMTSF)2PF6.^'^ In the phase diagram 
of (TMTSF)2PF6,^ the superconducting phase is adja- 
cent to the spin density wave (SDW) phase. From the 
knowledge of the anisotropic superconductivity,^"^ which 
have been extensively enriched for two decades, the d- 
wave singlet superconductivity (dSS) is expected to occur 
in such superconducting phase. However, recent experi- 
ments on (TMTSF)2PF6 indicate a counter example of 
this widespread knowledge. The unsaturated magnitude 
of the upper critical field Hc2^ and the constant Knight 
shift thorough the superconducting transition tempera- 
ture Tc^ strongly suggest the formation of triplet pairing. 
In addition, the SDW phase turned out to be an exotic 
one, in which the 2k p charge density wave (CDW) co- 
exists with the 2kF-SDW as seen from the X-ray exper- 
iments.^'* 

The theoretical study of qld systems is complicated 
due to its dimensionality originated from the anisotropy, 
where the effects of the one-dimensional (Id) fluctua- 
tion are essential to understand its physics. In the sys- 
tem with the Id fluctuation, there are equally diver- 
gent contributions from both Cooper (electron-electron) 
and Peierls (electron-hole) channels. However, the Id 
fluctuation cannot be evaluated correctly by the con- 
ventional method such as the random phase approxi- 
mation (RPA),^ the fluctuation exchange (FLEX) ap- 
proximation^^' or the third-order perturbation the- 
ory, which are powerful for the study of three- or two- 
dimensional anisotropic superconductivity. On the other 
hand, the renormalization group (RG) technique based 



on the g-ology is extremely useful for the treatment of 
the Id fluctuation, whereas it is complicated to deal with 
anisotropic superconducting states, such as dSS. 

Recently, Duprat and Bourbonnais extended the tech- 
nique of Kadanoff- Wilson RG for purely Id systems 
to that for qld systems. This technique enabled us 
to treat both Id fluctuations and various kinds of 
anisotropic superconductivity. Meanwhile, the essential 
role of charge fluctuation has been pointed out phe- 
nomenologically.^'^' Quite recently, it is shown by the 
RPA^^ that charge fluctuations induced by next-nearest 
repulsive interactions can give rise to the /-wave triplet 
superconductivity (/TS), though the Id fluctuation ef- 
fects are still pushed aside. 

The purpose of the present paper is to get a deep in- 
sight into the interplay of the superconductivity and the 
density waves in qld systems in terms of the Kadanoff- 
Wilson RG technique developed for qld system. Espe- 
cially, we shall focus on clarifying i) crossover from purely 
Id case into qld case, ii) effects of long-range Coulomb 
interactions and iii) microscopic mechanism of the qld 
superconductivity, considering the effect of the Id fluc- 
tuation. 

In §2, we briefly describe the Kadanoff- Wilson RG 
technique. Next, the method of the present work. In §3, 
it is shown that the response functions exhibit a no- 
ticeable property being proportional to cxp(l/r). Phase 
diagrams for several magnitudes of mterehain hopping 
are demonstrated on the plane of mtrachain interactions. 
Transition temperatures as a function of the mterehain 
hopping, are also presented. Then a discussion of the 
mechanism of the qld superconductivity is given. Sec- 
tion 4 is devoted to conclusion and a comparison of the 
present result and the experiments. 
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2. Formulation 

2.1 Model 

We consider a system consisting of an array of 
chains of length L, which is described by the Hamilto- 
nian, 



p,k,(7 
TTVF 



ki ,q,CT 
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= -TTVF W^ikib, k2b, Qb-, ^)0?*(q)0?(q), 

(10) 

where the renormalization of the bandwidth is defined as 

Eo{£)=EocM-i), (11) 

with the bandwidth cut-off Eq. The quantities ip(k) and 
tp*{k.) are the Grassmann fields and k = (k, a;„). {k}* 
denotes the renormalized momentum space. The opera- 
tor of. denotes Peierls fields 



^ "+,ki+q(Ti"-,k2-qCT2"-^k20-3a + ,ki<T4- (l) 

The first term denotes the band energy with the qld 
dispersion: 

^(k) = —2ta cos ka — 2tb cos kh — n 

- VF{\ka\ - kp) - 2tb COS kb + 2t't, COS 2kb, (2) 

where the inirachain electron spectrum with the nearest 
mirochain hopping ta is linearized at the Fermi points 
±kp with the Fermi velocity vp, and fx is the chemical 
potential. In Eq. (1), a^^^^ (opkcr) is a creation (an annihi- 
lation) operator close to the right ka = +kF {p = +) and 
the left ka = —kp {p = —) Fermi surface, respectively. 
The second and third terms of Eq. (2) denote the near- 
est mienchain hopping tb and the next-nearest mienchain 
hopping t[, respectively, where 

tl cos{k°p) 



t'u 



(3) 



4t„sin2(fc0)' 

is obtained as the nesting deviation. For the interaction 
part, the second term of Eq. (1), we follow the standard 
notations of the g-ology, where gi corresponds to back- 
ward scattering and (?2 does to forward scattering. (Note 
that gfi,2 are normalized by nvp here.) For the 1/4-filling 
case, gi and 92 are given in terms of on-site U , nearest- 
neighbor Vi and next-nearest-neighbor V2 interactions, 
as 



TTVpgi = U/2-V2, 
nvFg2 = U/2 + V1+V2. 



(4) 
(5) 



In order to examine superconductivity and density wave 

in terms of RG technique, we utilize the path integral 
method, where the partition function Z (= Tre"'^'^) is 
given by 

z = JJvrv^beMs[r,m, (6) 

and the effective action S is written in the form, 

s[r,^]i = So[r,^]i + Si[r,^]i, (7) 
so[r,^u= E [G°(k)]-v;.(k)Vp,.(k), (s) 

p,a,{k}' 

Si[r,i^]e = -nvF^J^ikn,k2b,qb;i)Ol*iQ)Ol{Q), 

(9) 



LNb 



E V*,a(k-Q)<Tf^^+,/3(k), 



«,/3,{k}* 



(12) 

describing the CDW = 0) and the SDW (/x = 1, 2, 3) 
channels, where Q = (2A:^ + Qq, Qb, ujm = 0, ±27rT, • • • ) 
(Qo < 2fc^). The operator O? is Cooper fields 



0^(q) 



LNb 



E ar,-a(-k + q)afV+,Mk), 



a,/3,{k}* 



(13) 



for singlet (/x = 0) and triplet (^ = 1,2,3) pairing, 
where q = (g'o,^,^^) (go < 2k^p). The couplings 
and Wfi correspond to that of the density wave and the 
superconductivity respectively, and are related to 51,2 as 
Jn=o = gi- 52/2, J;,/o = -52/2, W^=o = -(51 + 92)12 
and Wfi^o = {gi - g2)/2, respectively. 

2.2 Renormalization of couplings J^i.{(-), Wfi{(,) 

First we apply the RG technique to coupling constants 
for the Peierls and the Cooper channel. By taking con- 
traction of each channel, the one-loop flow equation are 
obtained as^^ 

^^Jn{Qb - kb,kb\t) 

= ^ E - , fc^; mnikl kb;£)lG{i) 



- MQb - kb, h; E - £)Ip{Qb, k'b, £), 
K 

(14) 

where Co,o = -l/2,co,^^o = 1/2,0^,^0,0 = 1/2 and 
c^^o.p^o = 1/6- The first term of the right-hand side 
of Eq. (14) corresponds to the contraction in the Cooper 
channel, where Ic{(-)d£ = tanh[i?o(^)/4T]d^ is the outer- 
shell integration. The second term corresponds to the 
contraction in the Peierls channel, with the outer-shell 
integration 

ip{Qb,kbi,i)dl = 



4 ^^Eoii)+XA{kbi,Qb) 



tanh 



Eoji) 
4T 



te.nh^{^+ XAikbuQb)} 
(15) 



A{kb, Qb) = 2i6[cos kb + cos(fcb -|- Qb)] 
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For tb = 0, the flow equation (14) yields the one-loop Id 
flow equations. -'^^ Note that /c(^) is a constant while 
IpiQb, kb, i) suppressed by i'j, depending on the wave vec- 
tor of particles. Thus the Cooper channel is always singu- 
lar (i.e., does not depend on the shape of the Fermi sur- 
face) while the Peierls channel is suppressed by the nest- 
ing deviation. The renormalized Cooper couplings are ob- 
tained from the renormalized Peierls couplings through 
the relation 



Wp,=Q{Qb -kb,kb]t) 



-i^JoiQb - kb, h'J) 
3 

-I- -Jn^o{Qb - kb,kb;i), 



Wfi^oiQb - kb,kb]l) = -JoiQb - kb,kb;i) 



+ -Jnjto{Qb - kb,kb'J), 



(17) 



(18) 



which is obtained from Eqs. (9) and (10) straightfor- 
wardly. 

Following the idea of the dimensional crossover the- 
ory/^' we adopt two-step RG procedure. First we apply 
the RG technique of Id^^ (Id RG) to the regime above 
the crossover temperature Taoss: and then that of qld^^ 
(qld RG) to the region below Tcross as the second step. 
For the first step, we adopt the two-loop fiow equations 
of ldi6 

dgi{e) 



dl 
dG{e) 



= 0, 



(19) 



(20) 



where G(£) = gi{e) - 2g2{t,. During the Id RG, the in- 
ienchain coupling j'^ due to the pair hopping is produced 
as follows: 



di 

dUM^) 

di 



ig2{e)-2g,{e))z{£)tb 



g2{i)z{i)tb 



n 2 



(21) 
(22) 

(23) 



where flow equations for tb and z (corresponding to self 
energy) are 

^ln4(^) = 1 -\[gi{e?+92{(?-9mgM^ (24) 



-.W?+g2{if-gi{i)g2{ 



(25) 



respectively. The dimensional crossover temperature is 
given by Tcross = i^o(4ross)/2, where 4ross is defined as 
-Eo(^cross) = 4(^cross)- Thus the initial couplings for the 
second step, i.e., qld RG, are given by a combination of 
5^ and cos Qb as 



J ii=o{,Qbi -Across ) 



-\[92{t. 



s) - 251 (4 



3 )] +ioCOS(56, 

(26) 



JtJ.jio{QbJcToss) = --^g2{icioss) + JnCOsQb- (27) 

2. 3 Response functions x 

In order to calculate the response function, we add the 
source field to Eq. (7), which is given by^^ 

Sh[r,i'] = ^z^[o;:*(Qp)/i^(Qp) +C.C] 
ft 

+ J2 4"^[o^*(k)/ilr'(fc6)+c.c], (28) 

Ai,{k}* 



where Qp = (2/c^,7r), and hp^\kb) = h''^^ cos{nkb) or 



hp^\kb) — hp'' sin(nfcf,). The corresponding pairing sym- 



,(0) 



metries are discussed later. Here ^^^"^^ corresponds to 
vertex corrections, and the value [2;^"!^]^ is equal to the 
well-known auxiliary susceptibility x-^^ Within the lin- 
ear response theory, only the linear and the quadratic 
term contribute to the action, so that the action at £ has 
the form 

V, h\ h]e = S[r,ij]e + J2 ^f^W{Ol*{Qp)h^{Qp) + cc] 

+ ^ 4"^(£)[Of (k)4")(fc6) + C.c] 

/i,{k}* 

H p,,ki, 

(29) 

Flow equations of each 2^"^) are given by 



di' 



lnz^{£) = -i-^ J^(7r- A;f„fc6;£)Jp(7r,fc6;^), (30) 



for the Peierls channel and 

±lnz^-2M = lcS^loii)Ici£), (31a) 

I ln45oW = ^cgoW^c(^), (31b) 

for the Cooper channel. Then the static response func- 
tions are obtained as 



The coefficients cS^ ^ represent the Fourier coefficients of 
the coupling in Cooper channel ajf^ and which are 
defined as 

Wfi{kb,k'b;i) =af\£) + J2[ap'\^)cos{nkb)cos{nk'b) 



n>0 



+ bf> {£) sm{nkb) sin(nA;;)] . (33) 

The coefiicicnts a^"' (f^"') are the even (odd) compo- 
nents with respect to the mierchain direction fc6,/c^. 
The parity of the gap cannot be classified only by kb- 
dcpendcncc. It has to be classified by the symmetry of 
both ka- and fcf,-dependences of the gap function. For sin- 
glet channel, considering only the even-frequency pairing, 
the parity of the gap is even, so that 
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Fig. 1. Types of the gap symmetry. The upper panels corresponds 
to the singlet pairing, whose symmetry belongs to s-, d^2_y2-, 
dxy-, g-wave respectively (from left to right). The lower panels 
corresponds to the triplet pairing, whose symmetry belongs to 
Px-, fx-, Py-, /y-wave respectively (from left to right). The co- 
efficients (a^"' , b^"' ) above the panels denote the corresponding 
Fourier coefficients. The arrows indicate the inversion (k — k), 
which does not (does) change the sign of the gap function for sin- 
glet (triplet) pairing. 



A) ap^2o' even-fca and even-fcf, (s, d^2_ 

Lin) 



•) 



B) odd-fca and odd-fc^ {dxy, g, ■■■) 

are possible. Similarly for triplet channel, the parity of 
the gap is odd, so that 

C) odd-ka and even-kb {px, fx 



■) 



D) b^pjo- even-/ca and odd-kb {py, fy, 
are possible. Briefly speaking, the symmetry of kb- 
dependence is classified by Op"'* or b^^\ and that of ka- 
dependence is automatically determined once we choose 
the spin-symmetry /2. 

In this paper, we consider the following 8 symmetries 
for order parameters of superconductivity. For singlet 
channel, we calculate the response functions for a^*Eo' 

flp^loi ^p=o ^"^^ ^;f=0' which can be recognized as s-, 
dx2-y2-, dxy- and g-wave, respectively.^^ For triplet chan- 
nel, the response functions for a^^^Q, ^^^^^o' ^fi^o ^^'^ ^/f#0' 
which can be recognized as p-, fx-, Py- and /^-wave re- 
spectively, are calculated. (See Fig. 1 for the graphical 
representation of each symmetry.) We examined whole 8 
symmetries, and found the clear instability only for px-, 
dx2_y2- and fx-wave. Therefore, from now on, wc simply 
call them p-, d-, /-wave, respectively. 

3. Superconductivity vs. Density Wave 

3.1 Temperature dependence of x{T) and Tc 

Based on the RG method derived in previous section, 
we examine the possible states in the present qld sys- 
tem. First we calculate the response functions obtained 
from Eq. (32), which are shown in Fig. 2 for the SDW 
(a) and the dSS (b). Each response function (also for 
other response functions) shows the similarity in behav- 
ior. The important point to note is that the response 
function x(T) becomes proportional to exp(l/T) (i.e., 



■) 




0.0 0.1 0.2 0.3 

i/r(K) 



0.4 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

i/r(K) 



Fig. 2. Response function of (a) SDW and (b) dSS. The dotted 
straight line is the line of Xprop drawn by fitting with x the 
inflection point. The insets show d{ln[iTVpx])/d{l/T) for the re- 
spective response function. 



Inx(r) cx l/T) in the wide region of intermediate tem- 
peratures. From the theorems of Mermin and Wagner,^" 
and Hohenberg,^^ there is no magnetic or supercon- 
ducting long-range order at any finite temperature in 
one and two dimensions. Generally, the response func- 
tions in two dimension have the temperature dependence 
x{T) oc exp(l/T) for both magnetic^'^ and superconduct- 
ing^-^"^^ cases, when we do not consider the Kosterlitz- 
Thouless transition. From these points, one may say 
that the present method gives a rigorous solution in such 
temperatures as x(T) cx exp(l/T) in spite of the one-loop 
approximation. 

In the very low temperatures, on the other hand, x(T) 
strays from the behavior of exp(l/T) and diverges at 
finite temperature because of the divergence of the cou- 
pling J or W. Therefore, it is appropriate to recognize 
the response function in the intermediate temperatures 
as a proper response function, since it is proportional 
to exp(l/T). In the region, where the obtained response 
function does not show the behavior of exp(l/r), the 
present RG technique does not work well due to the rel- 
evant (large enough) couplings. In the sense, the char- 
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acteristic temperature, T*, for the lower bound of the 
present RG technique, may be taken by a condition that 
d{hi[Trvpx\) /dT~^ takes a minimum, i.e., the inflection 
point. In order to obtain the proper response function, 
we estimate by using an extrapolation formula 



/,,= 100K 



(,,=2()()K 



Xprop(r) = aexp'^ 



(34) 



where x = ""^fX- The parameter /? is defined as /3 = 
d{ln[x]) / dT~^\T=T* , which is of the order of the transi- 
tion temperature obtained by the mean field approxima- 
tion. The parameter a is defined as a = x(T*)e~^^'^ , 
namely, Xprop(T') is the tangential line of lnx(l/r) at 
T = T*. Assuming that the practical transition tem- 
perature Tc is lead by the three-dimensionality, we use 
the conventional RPA (for the inter-plane coupling), 
53dXprop(T'c) = 1, which yields the form 

Tc = - , / (35) 

where g^d is an inter-plane coupling. Of course, we need 
to discuss the problem in the Kosterlitz-Thouless (KT) 
context in two-dimensional systems. In the case of 
Tc <C Ep, however, the difference between Tc and the 
KT transition temperature T^"^ is small compared to 
Te, i.e., |Te - TKT|/r, ~ Tc/Ep.^^ Consequently, de- 
termined by the present method almost corresponds to 
the actual transition temperature T^"^ . 

In this paper, Tc is calculated by taking 0{g3d) ^ 
10^^. We calculated 32 x 32 x 2 simultaneous differen- 
tial equations, i.e., a 32-chains system with 64 points 
on the Fermi surface. The initial bandwidth is set to be 
Eo = 5500K {ta = Eo/{2y/2k%) ~ 2500K ), which is 
consistent with the bandwidth of (TMTSF)2PF( 



28, 29 



3.2 g2-gi phase diagram 

First, we show phase diagrams on the g2-gi plane for 
several choices of t^'s in order to understand the variation 
of the phase diagram from pure Id case to qld case. For 
repulsive interactions 51,2 > 0, the SDW is dominant for 
gi < 2g2 and the p-wave triplet superconductivity (pTS) 
is dominant for gi > 2g2 as found in the g-ology."'^^ Figure 
3 shows the phase diagram for tb — 100, 200, 300, 400, 
500K. The dashed-lines in Fig. 3 denote 51=252, which 
correspond to the phase boundary between the pTS and 
the SDW of the g-ology. 

The present results show that the transition tempera- 
tures of the pTS and the SDW are extremely low (numer- 
ically, below 10~^ K) near the phase boundary gi — 2g2. 
Moving from that line, the pTS and the SDW phase ap- 
pears parallel to the gi=2g2-\ine. The salient feature is 
that, in qld, the dSS phase emerges between the gi=2g2- 
line and the SDW phase, and this dSS phase extends to 
the SDW phase also parallel to gi=2g2, as tb increases. 
From these features, it turns out that the transition tem- 
perature depends on the intensity of \gi — 2g2\, corre- 
sponding to the intensity of density fluctuations. The 
nearer the parameters move to the line of gi—2g2, the 
weaker the SDW instability becomes. Thus the dSS ap- 
pears instead of the SDW state due to the deviation of 
the nesting condition. Just on the line of 31=0, the CDW 
(/TS) state has the same transition temperature as the 











dSS 








SDW 



0.2 0.4 0.6 0.8 



Vo 

CDW. 
SDW 



pTS 



CDW, 
SDW 



/,,=300K 



/j=4(X)K 




Fig. 3. 92-31 phase diagrams for 4 = 100, 200, 300, 400 and 500K. 
The transition temperatures for the CDW, SDW, dSS, pTS and 
/TS are expressed as the tones of green, red, blue, yellow and or- 
ange, respectively. The dashed-line indicates the phase boundary 
between the pTS and the SDW of the g-ology. 



SDW (dSS) state. The phase boundary between the dSS 
and the normal state swerves from the line parallel to 
gi — 2(72 for small gi, while that between the dSS and 
the SDW is parallel to the gi—2gi-\me. This can be un- 
derstood as the interference effect of the charge fluctua- 
tions and the spin fluctuations, which shall be discussed 
later. 

3.3 tb-dependence of transition temperature 

Next, we show ifc-dependence of the transition tem- 
perature in Fig. 4; (a) the Hubbard {V=0) model with 
= 32 = 0.88 {U ~ 7.8ta,Vi = V2 = 0.0), and 
(b) the extended Hubbard model with large-y given 
by 51 = 0.01,52 = 0.55 {U ~ l.ita.Vi ~ 0.72^^ and 
V2 — 0.60ta, if we assume Vi = I.2V2). For small tb 
(<300K), the transition temperature of SDW, Tsdw, in- 
creases linearly. More precisely, it is approximately given 
by Tsdw ~ exp(— 2/52). Here, 52/^ corresponds to 
the strength of the spin fluctuation, where the asterisk 
denotes its value at Tcross- Such increase cannot be ob- 
tained by the mean field theory or the RPA-type the- 
ory, and consistent with the previous Id RG theory;^^''^" 
this feature originates from the effect of Id ffuctuation, 
i.e., the equally divergent contributions from both the 
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(a) 



g, = 0.88, g, = 0.88 
U=l\, 

yj=o.o\, V, = O.Or_^ 
£ =550dK,''« ~2500K 



SDW 
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(b) 



10 



/ \ 



\ 



=0.01, g, = 0.55 
C/= 1.3?^, 

yj=0.72r , K = 0.60r 
£ =55Qok./~2500k 
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(CDW) 



(/rs) 
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400 600 
^.(K) 



800 



Fig. 4. Transition temperature as a function of the mterehain 
hopping t\, for (a) the Hubbard (y=0) model gi = 52 = 0.88 
{V ~ 7.8ta,Vi = V2 = 0.0), and for (b) the large-V model 
51 = 0.01,52 = 0.55 (U ~ 1.3ta,Vi ~ 0.72ta,V2 ~ 0.60ta). 
The closed circles and triangles are the transition temperature 
of SDW and dSS, respectively. The open circles and diamonds 
indicates the temperature where the response functions of CDW 
and /TS satisfy the condition (35), respectively. 



Cooper and the Peierls channel. The transition tempera- 
ture has a maximum followed by the decrease due to the 
nesting deviation. When TgDw is sufficiently suppressed 
due to the nesting deviation, the hidden superconductiv- 
ity, dSS in the present case, emerges. These features are 
seen for almost every couplings in the parameter region 
< gi < 2g2, corresponding to the SDW phase of the 
g-ology. 

There are some differences between the V—Q and the 
large-F model. The most noticeable difference exists in 
the subdominant phase for the large-1^, which is shown 
by open squares (CDW) and diamonds (/TS) in Fig. 4 
(b). Here, the transition temperature for the subdom- 
inant state are also obtained from the condition (35). 
The transition temperature for the subdominant CDW 
(/TS) is exactly the same as that of SDW (dSS) when 
(?i = (f/ = 2V2), i.e., just on the phase boundary be- 
tween SDW and CDW of the g-ology. The actual phase 



transition, of course, occurs only for the most divergent 
response function, so that these subdominant phases are 
masked practically. Therefore the subdominant transi- 
tion temperature is better considered as an emergence of 
fluctuations. In the case of large- 1^, not only the spin fluc- 
tuations but also the charge fluctuations develop, while 
only the spin fluctuations develop for the V=0. This en- 
hancement of the charge fluctuation promotes the /TS 
state; especially, /TS is realized for gi — 0. 

Another significant difference is that the supercon- 
ducting transition temperature Tc for large-y is much 
lower than that for 1^=0, while Tsdw of both cases are 
almost the same. This comes from the difference of the 
strength of the charge fluctuations, ft is clear from the 
eqs. (17) and (18) that the charge fluctuation Jq sup- 
presses the singlet channel Wq, and enhances the triplet 
channel Wp^^Q. When Jg = Jfj,^o, the amplitude of the 
triplet channel becomes just the same as that of singlet 
one, which corresponds to the 51 = (C/ ~ 2V2) case. 
Thus the low Tc of the large- model also suggests the 
enhancement of the charge fluctuation, which promotes 
the /TS. The tb-dependence of Tc would be given by 
Tc ~ tbexp(— l/g^) in the same way as the SDW. (Here, 
g*^ is slightly smaller than g*^ because of the nodes of the 
gap.) Therefore, the increasing Tc for the Hubbard model 
is due to the increasing if,. The reason why Tc decreases 
for the large- is the rapid decrease of the spin fluctua- 
tion due to the cancellation with the charge fluctuation. 

Note that the possible /TS state is a characteristic 
of systems with open-Fermi-surface such as qld. On the 
other hand, in the case of the closed Fermi surface, such 
as the square lattice, /TS cannot have the compara- 
ble transition temperature even for large-T^, but d^y- 
wave singlet superconductivity can become a dominant 
state. '^'^''^^ This comes from a property of the open Fermi 
surface. The number of nodes of /TS is 4 in the case of 
qld open Fermi surface, which is the same as dSS. Here, 
the number is 6 for the closed Fermi surface. Hence, the 
qld system possesses a property which is favorable for 
the triplet pairing. As far as the pairing symmetry is 
concerned, the present results is consistent with the RPA 
results. The RPA, however, does not take into account 
the Id fluctuations, so that its Tc is much higher than 
the present one. 

4. Conclusion 

We have calculated the response functions and the 
transition temperatures for the density wave (CDW 
and SDW) and the superconductivity (s-, Px-, Py-, 
dx2-y2-, dxy-, fx-, fy- and g-wave) of qld systems by 
the Kadanoff- Wilson renormalization group technique, 
which has been developed by Duprat and Bourbon- 
nais.^^ The present method is able to investigate a sys- 
tem in pure Id and qld in the same footing. Each re- 
sponse function exhibits the temperature dependence, 
x(T) oc exp(l/T), in a wide region of temperature, which 
satisfies the theorems of the Hohenberg,^^ and Mermin 
and Wagner. ^° Also, the temperature dependence is con- 
sistent with that obtained by 2d fluctuation theories. ^^"^^ 
This may suggest that the present method gives a rigor- 
ous solution to some extent. 
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In the phase diagram of the (72-ffi plane, the pTS and 
the SDW are locating in the region parallel to gi=2g2. 
The transition temperature depends on the intensity 
Igi — 2(72 1, which corresponds to the coupling constant 
of density- fluctuations. The dSS phase, which does not 
appear in purely Id case, exists next to the SDW phase 
in qld case. The dSS phase expands into the SDW re- 
gion parallel to g\=2g2 as t;, increases. For small 
TsDw increases as Tgovv ^ cxp(— 2/52)- The transi- 
tion temperature of the SDW has a maximum and de- 
creases at finite % (~ 300K for (51,52) = (0.88,0.88) 
and (51,52) = (0.01,0.55)), due to the nesting devia- 
tion. The singlet superconductivity of (ix2-j/2-wave (dSS) 
emerges when the SDW state is suppressed by the mod- 
erately largo tb. In the large-F model, the subdominant 
CDW (/-wave triplet) state is comparable to the dom- 
inant SDW (da,2_y2-wave singlet) state. Especially, the 
transition temperature of /TS (CDW) becomes equiva- 
lent to that of dSS (SDW) with V2 = U/2. This is because 
the large V (practically V2) leads to the enhancement 
of the charge fluctuations. Such large charge-fluctuation 
suppresses the dSS and enhances the /TS, so that 
of dSS in the large-V^ model is much lower than that in 
the V=Q model. Consequently, the long-range Coulomb 
interaction can give rise to the large charge fluctuation, 
which suppresses dSS and enhances /TS particularly in 
qld systems. 

The experiments on (TMTSF)2PF6 compounds sug- 
gest the coexistence of 2A:i7-CDW and 2fci?-SDW,^'^ and 
exhibit about ten times lower Tc than Tsdw-^ Consid- 
ering these experimental properties, the \aige-V model 
{y2/U ~ 0.5) is appropriate for the (TMTSF)2PF6 com- 
pounds on the point of the phase diagram, the 2kF- 
charge-fluctuation in SDW phase, and the realization of 
the /TS. Actually, there is an estimation that V\/U ~ 
0.5-0.75 {V2/U ~ 0.4-0.6, if we assume Vi = I.2V2) 
on a basis of a Valence-Bond/Hartree-Fock method for 
molecules derived of TTF.^^ 

Even if V2 is not so largo, /TS is possibly realized 
under moderately large Zeeman magnetic field, which 
suppresses the singlet pairing of dSS due to the param- 
agnetic effect. In such a case, since the ground state 
is (iSS state, the Knight shift decreases at H = 
but does not vary under large magnetic field. Noting 
that the upper critical field can be estimated as Hc2 = 
(Ao/\/2/iB)(l - TVr<?),34 we can evaluate that 1.7 T is 
sufficient to realize the field-induced /-wave triplet pair- 
ing for Tc{H = 0) = 1.2K, even for T/ ~ T;^/2.^^ Here, 
Aq is the magnitude of the superconducting gap, and the 
BCS relation 2Ao = 3.5A;bT'c is assumed. 

Finally, we comment on the effects of the dimerization 
and the Umklapp scattering, which are not considered 
in the present work. Following effect is expected with 
the Umklapp scattering 53 as tb increases. ^^'^^ First the 
Umklapp scattering causes the antifcrromagnctic transi- 
tion (53 is relevant while 51 is irrelevant). Next, the SDW 
transition is realized due to the nesting property (51 be- 
comes relevant). Finally, the /TS state results from the 
charge fiuctuation enhanced by the off-site interaction, 
which is not affected by the Umklapp scattering; the es- 
sential points of the present paper would be still valid. 
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